Multigrid with rough coefficients and Multiresolution operator decomposition from Hierarchical Information Games
نویسنده
چکیده
We introduce a near-linear complexity (geometric and meshless/algebraic) multigrid/multiresolution method for PDEs with rough (L∞) coefficients with rigorous a priori accuracy and performance estimates. The method is discovered through a decision/game theory formulation of the problems of (1) identifying restriction and interpolation operators, (2) recovering a signal from incomplete measurements based on norm constraints on its image under a linear operator, and (3) gambling on the value of the solution of the PDE based on a hierarchy of nested measurements of its solution or source term. The resulting elementary gambles form a hierarchy of (deterministic) basis functions of H1 0 (Ω) (gamblets) that (1) are orthogonal across subscales/subbands with respect to the scalar product induced by the energy norm of the PDE, (2) enable sparse compression of the solution space in H1 0 (Ω), and (3) induce an orthogonal multiresolution operator decomposition. The operating diagram of the multigrid method is that of an inverted pyramid in which gamblets are computed locally (by virtue of their exponential decay) and hierarchically (from fine to coarse scales) and the PDE is decomposed into a hierarchy of independent linear systems with uniformly bounded condition numbers. The resulting algorithm is parallelizable both in space (via localization) and in bandwidth/subscale (subscales can be computed independently from each other). Although the method is deterministic, it has a natural Bayesian interpretation under the measure of probability emerging (as a mixed strategy) from the information game formulation, and multiresolution approximations form a martingale with respect to the filtration induced by the hierarchy of nested measurements.
منابع مشابه
On Multiresolution Methods in Numerical Analysis 4832
As a way to emphasize several distinct features of the mul-tiresolution methods based on wavelets, we describe connections between the multiresolution LU decomposition, multigrid and multiresolution re-duction/homogenization for self-adjoint, strictly elliptic operators. We point out that the multiresolution LU decomposition resembles a direct multigrid method (without W-cycles) and that the al...
متن کاملA Multigrid Approach to Image Processing
A second order partial differential operator is applied to an image function. By using a multigrid operator known from the so-called approximation property, we derive a new type of multiresolution decomposition of the image. As an example, the Poisson case is treated in-depth. Using the new transform we devise an algorithm for image fusion. The actual recombination is performed on the image fun...
متن کاملLU Factorization of Non-standard Forms and Direct Multiresolution Solvers
In this paper we introduce the multiresolution LU factorization of non-standard forms (NS-forms) and develop fast direct multiresolution methods for solving systems of linear algebraic equations arising in elliptic problems. The NS-form has been shown to provide a sparse representation for a wide class of operators, including those arising in strictly elliptic problems. For example, Green’s fun...
متن کاملMultiresolution Approximate Inverse Preconditioners
We introduce a new preconditioner for elliptic PDE’s on unstructured meshes. Using a wavelet-inspired basis we compress the inverse of the matrix, allowing an effective sparse approximate inverse by solving the sparsity vs. accuracy conflict. The key issue in this compression is to use second-generation wavelets which can be adapted to the unstructured mesh, the true boundary conditions, and ev...
متن کاملFiedler trees for multiscale surface analysis
In this work we introduce a new hierarchical surface decomposition method for multiscale analysis of surface meshes. In contrast to other multiresolution methods, our approach relies on spectral properties of the surface to build a binary hierarchical decomposition. Namely, we utilize the first nontrivial eigenfunction of the Laplace-Beltrami operator to recursively decompose the surface. For t...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM Review
دوره 59 شماره
صفحات -
تاریخ انتشار 2017